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We study the possibility of performing quantum state tomography via equidistant states. This 
class of states allows us to propose a non-symmetric informationally complete POVM based tomo- 
graphic scheme. The scheme is defined for odd dimensions and involves an inversion which can be 
analytically carried out by Fourier transform. 



I. INTRODUCTION 

The determination of quantum states and processes 
[l| plays an important role in the foundation of quantum 
theory 0, @ and quantum information theory J5,, J,, ^8,, . 
Quantum states are mathematical descriptions of physi- 
cal systems, which provide all the information necessary 
for characterizing the result of any possible experiment. 
According to this, quantum states, unlike their classi- 
cal counterparts, are not observable; therefore the de- 
termination of quantum states becomes an important 
problem. Furthermore, in quantum information theory 
the implementation of processes requires a quality as- 
sessment, which typically bases on the determination of 
quantum states. 

The most commonly used method to determine 
unknown quantum states belonging to a finite N- 
dimensional Hilbert space is quantum tomography [Tol . 
[llj . This method requires the preparation of a large en- 
semble of copies of the unknown state and the capacity of 
measuring transition probabilities toward a set of states, 
which when properly chosen leads to the determination of 
the unknown quantum state. In general, the coefficients 
of the unknown state are lineal functions of the transition 
probabilities and thus an inversion is required. Due to 
this inversion and the inherent experimental inaccuracies, 
the reconstructed mathematical object might violate the 
positivity constraint defining quantum state. To over- 
come this problem, techniques like maximum likelihood 
estimation 12, 13'] or Bayesian analysis are applied. 

Several schemes for quantum tomography have been 
proposed depending on the choice of the set of states. 
Standard quantum tomography bases on the measure- 
ment of transition probabilities, calculated according to 
Born statistical formula, toward the N eigenstates of a 
A^-dimensional representation of the N'^ — 1 generators 
of SU{N) 14|, |l5|. A second alternative is symmetric 
informationally complete (SIC) POVM based quantum 
tomography [11, [13, [H, which resorts to the mea- 
surement of transition probabilities toward a set of N"^ 
linearly dependent states such that the absolute value 
of the inner product between two of them is 1/ \/N + 1. 
The existence of SIC-POVM has been shown analyti- 
cally for iV < 8 and numerically for < 47. The ex- 
istence of SIC-POVM in arbitrary, finite dimensions has 
been conjectured but not proven. A third alternative is 



quantum tomo grap hy based on mutually unbiased bases 
(MUB) [20l. [2ll.l22j. which resorts to the measurement of 
transition probabilities to states belonging to A^ -I- 1 A^- 
dimensional orthonormal bases. States belonging to dif- 
ferent bases have an absolute value of the inner product 
equal to 1///V [13, Hi]. It has been shown that MUBs 
exist when N is an integer power of a prime. Otherwise 
the existence of MUBs is an open problem. 

In this article we introduce a new scheme for imple- 
menting quantum tomography. This is based on the so- 
called equidistant states [251]. These are such that the 
inner product between two of them is a complex con- 
stant or its conjugate. Projectors onto these states allow 
us to define a complete POVM which reconstructs an 
unknown density operator acting on a Hilbert space of 
finite, odd dimension A. The POVM can be analyti- 
cally constructed for these dimensions being the number 
of measurements to be performed equal to A^, as few as 
in the case of SIC-POVM or MUB-POVM based quan- 
tum tomography. The scheme requires inversion. How- 
ever, this is reduced to the inversion of (A^ -I- l)/2 circu- 
lant N X N matrices, which can be performed by Fourier 
transform. Thereby, this new scheme lays between stan- 
dard quantum tomography and POVM based quantum 
tomography. 



II. EQUIDISTANT STATES 



Let us start by reviewing the set Bo{ct) of equidistant 
states. The A non-orthogonal states \aj) {j = 0, . . . , N — 
1) in this set belong to a A^-dimensional Hilbert space Ti 
and are defined by the property 



(1) 



that is the inner product between any two equidistant 
states is equal to a or its conjugate a*. The states in 
Bo{a) are linearly independent when the modulo of a 
is constrained to the interval [0, |ae|), where \ae\ is a 
function of 9 and A^ given by 
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which is the largest possible value of the inner product 
|a| for and A^ fix. The states in yBo(a) become linearly 
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dependent when \a\ = \ag\ for 9^0, lying symmetrically 
on a {N — l)-dimensional subspace. In the case 6 = 
the value of the bound \ag\ becomes one. Consequently, 
states with 9 = and are all within the same ray. 

The set of equidistant states So (a) becomes a set of 
symmetric states when the inner product a is real [25j . 
that is 6* = or 6* = TT. For these values of 9 the upper 
bound \ag\ we obtain jaol = 1 and jajrl — ^/{N — 1) 
respectively. In general, for a complex inner product a 
equidistant states are different from symmetric states. 

A canonical decomposition of equidistant states [26j | is 
the following 



(3) 



The coefficients Afc are eigenvalues of the matrix M con- 
taining the inner products {ai\aj). These are given by 



Afc = 1- lal 



N 
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and fulfill the identity 



and thus the iV^ projectors 11* form a positive opera- 
tor valued measurement. Thereby, the set of projectors 
seems to be well suited for reconstructing an unknown 
density matrix. 

We now consider an arbitrary unknown density oper- 
ator p acting onto Ti. This space is spanned via the 
canonical base {|fc)} with k = . . . N — 1. Thereby, the 
density operator can be represented as 



P = 



= J2 Pk,k'\k){k'\ 



(10) 



fc,fc'=0 



The probability Pj" of projecting the density operator p 
onto the state j-th state jap of the base Bs{a) is calcu- 
lated through the Born statistical formula 



According to Eq. ^ the set of N'^ probabilities Pj obeys 
the constraint 
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The complex phases ujk are the coefficients entering in 
the matrix T which diagonalizes the matrix M. These 
phases are given by 
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III. TOMOGRAPHY 



In order to reconstruct the state of the system we con- 
sider the sets Bs {a) with s = ... TV — 1 defined as 



(7) 



where the action of the unitary operator X onto the 
canonical base is 



X\k) = 1), 
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and the addition is performed modulo N . Let us note 
that X^ = X^ = 1 and that the unitarity of X indicates 
that each set Bs{a) is composed of equidistant states. 
Furthermore, all sets Bs{ol) are characterized by the same 
value a of the inner product. 

Now, defining the projectors IK = |a|)(aj| onto 
equidistant states and using Eqs. ([Hi, ([5]) and it can 
be shown that 
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Considering the decomposition of the equidistant states 
onto the canonical base the probabilities P/ become 
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(13) 



p,g=0 



where the operations with subindexes are carried out 
modulo N . 

The previous expression leads us to a linear equations 
system with N"^ unknown variables, the real and imagi- 
nary parts of the coefficients of the density operator p 
spanned on the canonical base. Since the number of 
equations, which is given by the number of measurable 
probabilities PJ , is exactly iV^ {N bases x N states) the 
system can be inverted, provided the determinant of the 
transformation does not vanish, and consequently, the 
density operator can be reconstructed. 

However, the case N even presents a peculiarity. Cast- 
ing Eq. in the form 



N 
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where the coefficients pp^q have been separated in its real 
^{Pp,q) and imaginary 3(pp g) parts, we note that the 
coefficient {p — q) entering in Eq. (fTi|) assumes integer 
values between N —1 and 1, in particular the value N/2. 
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Thereby, the argument entering in the sine function in 
Eq. becomes an integer multiple of tt for any value 
of j and s. Thus, the imaginary part of the coefficients 
Pp,q with {p — q) = N /2 does not appear in the equations 
system and the previous set of probabilities does not al- 
low us to reconstruct the state completely. Therefore, 
the previous sets of equidistant states are at most well 
suited for reconstructing density operators acting onto 
Hilbert spaces with N odd. In the following sections we 
show that this is the case. 



IV. CASE N=3 

Let us now study Eq. ([TH)) in the particular case A'' = 3. 
In order to cast this equation in a simpler form we resort 
to the Fourier transform. Thereby we obtain 



where 



N-l 
q=0 
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Considering the three possible values of s we generate the 
following equations system 
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This can be separated into two independent sets of equa- 
tions considering the real and imaginary parts of the co- 
efhcients pp^q and of Pf. These two equations systems 
are characterized by the same matrix, which in our case 
turns out to be a circulant matrix, that is each row of 
the matrix is a cyclic shift of the previous row, being the 
eigenvalues given by the Fourier transform of one of the 
rows. Thus, the equations systems can be analytically 
inverted. 

The diagonal coefficients pk,k can be obtained by not- 
ing that 
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where 



(18) 
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which turns out to be also defined by a circulant matrix. 

Equations systems of Eqs. pT)) and (PO)) can be con- 
siderably simplified by taking Ai = 0. This particular 
choice is equivalent to consider the equidistant states in 
the set Bo{a) to be linearly dependent. In this case we 
have simpler solutions for the non-diagonal coefRcients 
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and for the diagonal coefficients we obtain 
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The choice Ai = is equivalent to make \a\ = jaej. 
It is still possible to select a particular value of \ag\ by 
choosing the value of the phase 9. In the particular case 
6 — TT the inner product a becomes —1/2, Aq = A2 = 3/2, 
(luo)* = e"^'''/^ and (^2)* = e^'^'^^. Thereby, the states 
in the sets Bo{a), Bi{a) and 82(0) are correspondingly 
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These states are such that the modulo squared of the 
inner product between any two of them is 1/4 and con- 
sequently they form a SIC-POVM. 



V. CASE N ODD 



The three possible values of s lead to the equations sys- 
tem 
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The previous results lead us to consider the following 
transformation of the transition probabilities 
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which when apphed onto Eq. ((T3)) lead us to the following 
equation 



N-l 
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(25) 
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where the eigenvalues 7^'^'' are given by 
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Tr 
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For fc and s fixed, this equation contains exactly dif- 
ferent coefficients of the density operator. The N values 
ofs = 0,...,iV— 1 lead to a linear equations system, 
which for a given k involve the N coefficients of the k- 
th diagonal of p. In particular, the equations system for 
k = involves the coefficients of the main diagonal only. 
The matrix M^*^^ defining the equations system has real 
coefficients only, thus it is possible to generate two equa- 
tions systems each for the real and imaginary parts of 
the coefficients of the fc-th diagonal of p. Thereby, we 
obtain one equations system for the N real coefficients of 
the main diagonal of p and {N — l)/2 equations systems, 
each for N complex coefficients belonging to a particular 
diagonal of p. This allow us to obtain the value of the iV^ 
real coefficients determining the unknown density matrix 
P- 

The matrix Af C') is circulant for A: = 0, . . . , (iV - l)/2, 
that is every row of the matrix M^*^' is a right cyclic shift 
of the row above. This class of matrices can be expressed 
in the form M^'^^ — F^D^'^^F, where F is the discrete 
Fourier transform and Z?'''"' — diag{'^'l^\^^^\ •■■i7^1i) is 
a diagonal matrix which contains the eigenvalues of M^'^^ . 

(k) 

The eigenvalues 7^ can be written as functions of the 
coefficients M^^} of the matrix M^^', with t = 0, iV- 1 
in the following way 
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The inverse of the matrix M^''\ provided M^*^^ is non- 
singular, is also circulant and is given by [mW]-i = 
i^t[i:)(fc)]-ii^^ where in this case the inverse of the dC') 
matrix is given by 
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Therefore, the elements of the first row that defines 
are given by 



fc)i 



. N-l , (fe)w 
r=0 I7i- I 



0,...,iV-l. (28) 



Let us consider for instance the case = 0. In this case 
we obtain that the values of are given by 



Ao = l + (A^-l)|a| 
Afc = fc=l, 

(fe) 
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and thus the eigenvalues 7^ vanish only for \a 
that is when the equidistant states in Bq(q) are mutually 
orthogonal. A similar result holds in the case 6 — tt. 



VI. CONCLUSIONS 

In the previous sections we studied the possibility of 
applying equidistant states to the problem of determining 
unknown states of finite dimensional quantum systems. 
We have shown that when the dimension N of the quan- 
tum system is odd, it is possible to formulate a scheme 
for quantum tomography. This scheme has two attributes 
that place it between standard quantum tomography and 
POVM-based quantum tomography. Firstly, the new 
scheme requires the measurement of iV^ transition prob- 
abilities, as few as in the case of MUB or SIC quantum 
tomography. Secondly, the new scheme requires the solu- 
tion of a system of linear equations. Thus, it requires to 
perform an inversion as in the case of standard quantum 
tomography. However, this inversion reduces to the prob- 
lem of inverting {N + l)/2 matrices each one of iV x TV. 
Since these matrices are circulant analytical expressions 
are known for their eigenvalues and eigenvectors. 

An interesting feature of the equidistant states based 
quantum tomography is that each one of the (TV -|- l)/2 
matrices provides information about a particular diag- 
onal of the density matrix to be determined. Thereby, 
errors in the reconstruction are constrained within diag- 
onals and do not propagate to other diagonals. Finally, 
let us note that the value of a can be freely chosen as 

(k) 

long as the eigenvalues 7^ do not vanish. 



Considering Eq. (|25p we obtain that the coefficients of 
the density matrix are related to the probabilities Pf 
through the relation 



1 



Pk+q,q 



N-l 
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TV ^ uWp 



'-[{l-q}r+kj] pi 



(29) 
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